A New Weak Lensing Analysis of MS1224.7+2007 



Philippe Fischeif] 
Dept. of Astronomy, University of Michigan, Ann Arbor, MI 48109 

ABSTRACT 

"on : 

ON ■ 

On ; Galaxy cluster mass distributions are useful probes of fl and the nature of 

the dark matter. Large clusters will distort the observed shapes of background 
galaxies through gravitational lensing allowing the measurement of the cluster 
mass distributions. For most cases, the agreement between weak lensing and 

C\| \ radial velocity mass measurements of clusters is reasonably good. There is, how- 

ever, one significant exception, the z = 0.32 cluster MS1224+2007, which has a 

' lensing mass substantially larger than the virial mass and also a very high mass- 

j — ' 

to-light ratio. Since this controversial object might be an unusually dark mass a 
follow-up study is definitely warranted. 
<^> In this paper we study the mass and light distributions of MS 1224+2007 out 

to a projected radius of 800 h^ 1 kpc by measuring the gravitationally-induced 
' distortions of background galaxies. We detect a shear signal in the background 

Q-f galaxies in the radial range 27.5" < r < 275" at the 5.5a level. The resultant 

mass map exhibits a peak centered on the dominant cluster galaxy and strong 

^ ' evidence for substructure which is even more strongly seen in the galaxy distribu- 

tion. Assuming all the detected shear is due to mass at z = 0.32 we find cluster 
mass-to-light ratio of M/L^ = 640 ± 150 (M/Lr) . The mass profile is quite flat 
compared to other clusters, disagreeing with a pseudo-singular isothermal sphere 
at the 95% confidence level. Our mass and M/L estimates are consistent with 
the previous weak lensing result. The discrepancy between the lensing and virial 
mass remains although it might be partially explained by subclustering and infall 
perpendicular to the line-of-site. This cluster remains a candidate dark object 
deficient in baryons and as such severely tests cosmological models. 

Subject headings: dark matter — galaxies: clusters: individual (MS1224. 7+2007) 
— gravitational lensing 



1 Bascd in part on research carried out at the MDM Observatory, operated by Columbia University 
Dartmouth College, University of Michigan and Ohio State University 
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1. Introduction 

Weak lensing distortions of background galaxies are a powerful tool for mea- 
suring galaxy cluster mass distributions. Lensing based studies are direct mea- 
surements and are not model-dependent as are other techniques (X-ray, radial 
velocities). Currently there are around 20 clusters which are well-studied with 
lensing (e.g. |Fischer fc Tyson (1997)] ), most yielding an evolution corrected 



M/L« 300 (M/L) . 

There is, however, one significant exception, the z=0.32 cluster MS1224+2007, 



which has M/L = 800 (M/L) from the weak lensing measurements of Fahlman 
et al. (1994)| . Not only is the lensing- derived M/L estimate very high but the 
mass estimate is approximately a factor of 2.5 higher than the virial mass esti- 
mate flCarlbcrg et al. 1994)) . If the lensing mass is correct then it would indicate 
a large region with an anomolously low baryon fraction. Within the standard 
cosmological model there is no causal mechanism which can generate primordial 
fluctuations in the baryon-to-total mass ratio on cluster scales ( Lvrard 1997| ). 



Furthermore, dynamical processes operating differentially on the baryonic and 
dark matter do not appear able to cause variations at the implied level (i.e. 
Metzler fc Evrard 1994]) . 



In general, clusters are found based on their optical appearance and/or their 
X-ray emission, which are both baryonic in origin. However, baryons are a small 
component of the total cluster mass (0.060±0.003/i- 3 / 2 , |Evrard 1997|) . Therefore, 



clusters may be very biased locations; measurements made in these regions may 
not be representative of the whole universe. Are clusters special places where 
there just happens to be large amounts of baryonic matter? Are there massive 
objects containing little or no baryonic matter which we have yet to detect? Is 
MS1224 an example of such a dark object? If so, this would represent a significant 
challenge for theories of large-scale structure formation and might also imply that 
estimates of Q based on cluster masses are underestimated. 

It is, therefore, important to confirm the weak lensing mass measurement 
of MS 1224+2007 with an independent study. In this paper we present a weak 
lensing study of MS 1224+2007 carried out with the MDM 2.4m. The observations 
and reductions are described in §[2] and §[| The shear and mass measurements are 
described in §[| and §[| along with comparisons with previous lensing and virial 
mass estimates. Conclusions and future work are presented in §| 
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2. Observations 

The cluster MS 1224+2007 was observed using with Michigan-Dartmouth- 
MIT (MDM) 2.4m telescope on 8-9 Feb. 1998. The total exposure time was 
21600s in R (41 x 900s) and 4500s in B (5 x 900s). The telescope was dithered 
between exposures. The "Echelle" 2048 2 thinned SITe CCD was used with 0.275" 
pixels. The seeing on the combined R image is 0.9" FWHM and 1.4" for the com- 
bined B image. The RMS sky noise values are 28.0 B mag per square arcsecond 
and 28.1 R mag per square arcsecond. The two nights were not perfectly photo- 
metric, the sky transparency varied by 9% as determined from monitoring 50 stars 
in the field of MS 1224+2007 (due to occasional small clouds passing through the 
field of view). We scale each image to match the flux of the highest throughput 
images. 

Despite the variations in sky transparency, 50 standard stars flLandolt 1992| ) 
were measured over the two nights. Transformation functions with linear color 
and airmass terms are fit to the data resulting in an RMS of 0.03 mag and 0.013 
mag for the B and R data respectively while the color terms are 0.035 and 0.013, 
respectively. Since the R-band data is much deeper than the B-band data and 
the color terms are very small, we calibrate the data assuming a mean galaxy 
color of (B-R)=1.6. This introduces errors of ±0.05 and ±0.02 for the B and R 
photometry respectively. 

The reddening in this field is E(B-V)=0.04 yielding A# = 0.17 and Ar = 0.11 
( [Schlegel et al. 1998| ). 

3. Faint Galaxy Photometry and Analysis 

The faint galaxy analysis was carried out using the analysis software ProFit 
(developed by the author). This software, starting with the brightest detections, 
fits an analytical model to each object using weighted, non-linear least squares, 
and subtracts the light from the image. It then proceeds to successively fainter 
objects. Once it has detected and subtracted all the objects in an image it 
replaces each in turn and refits and resubtracts until convergence is achieved. The 
software outputs brightness, orientation, ellipticity and other image parameters 
based on the fitted function. 



4. Gravitational Shear 
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4.1. Theory 

For gravitational lensing, the relationship between the tangential shear, jt, 
and surface mass density, E, is ( [Miralda-Kscude 1991| , |Miralda-Kscude 1995| ), 



lT\r) = «(< r) - «(r), (1) 

where k = S/E cr j 4 , the ratio of the surface density to the critical surface density 
for multiple lensing, and r is the angular distance from a given point in the 
mass distribution. The critical density depends on the redshift distribution of 
the background galaxies. The first term on the right is the mean density interior 
to r and the second term is the mean density at r. Therefore, the presence of a 
foreground mass distribution will distort the appearance of background galaxies. 
For a given coordinate on the image (f) , the distortion quantity for the i th galaxy 
is: 



r>^ _ 1-Oi/Qi) 2 [cos(2^)(Ax, 2 - A Vl 2 ) +2sm(26 i )Ax i Ay i ] 
l[r) ~ 1 + {k/atf X Axf + Ay* ' { ) 

where (fej/cij) and 6i are the galaxy axis ratio and position angle, respectively. 
Ax and Ay are the horizontal and vertical angular distances from f to galaxy i. 
D is related to the tangential shear by flSeitz fc Schneider 1995Q : 

< D(r) >= 2 m 1 \ n 3 

In the weak lensing regime k << 1, and 7r << 1, 7t ~< -D > /2. 

Before proceeding to a discussion of the shear measurments we discuss sources 
of systematic errors. 



4.2. Point-Spread Function Anisotropy 

In Figure [I] the ellipticity and orientation for 87 stars with 16.4 < R < 22.3 in 
the combined R-band image are shown. There is a strong anisotropy present in 
the point spread function (PSF) which, if not corrected for, could bias estimates of 
the shear due to the cluster. We correct for the PSF anisotropy using the method 
outlined in [Fischer fc Tyson (1997)| , which involves deriving a position dependent 
kernel which, after convolution with the image, yields round PSFs. Figure [l] 
shows the ellipticity and orientation for the same stars after the correction. We 
discuss the effects of PSF anisotropy on the shear and cluster mass estimates in 
50 and 953 
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Fig. 1. — Ellipticity and orientation for 87 stars in the combined R-band image. The lower 
panel is the original image and the upper panel is after correction (see text). The maximum 
ellipticity is e = 0.07 and the medians are < e >= 0.035 and 0.013, for the original and 
corrected, respectively. North is up and East is to the left. 
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4.3. Seeing and Shear Polarizability 

In order to measure the cluster mass using weak gravitational lensing we 
must measure the shapes of background galaxies to estimate the shear (see § |4.1|) . 
Even perfect galaxy measurements will result in biased shear estimates for two 
reasons. The first is that the galaxy ellipticities will be underestimated due to 
blurring by atmospheric seeing. The second is that the response of a galaxy to 
a shear will depend on its ellipticity and orientation with respect to the shear 
and hence the mean response will depend on the intrinsic ellipticity distribution 
of the galaxies. In order to calibrate these effects we carry out simulations using 
the F606W Hubble Deep Field Data (HDF) QWilliams et al. 1996Q , using the 
techniques described in [Kaiser et al. ( 1995 )| . This involves stretching the HDF 
data by 1 + 5, convolving with the PSF and adding noise. The values of Di (see 
Equation 0) are measured for each galaxy and compared to the unsheared values 
of Di. The quantity of interest is the recovery factor, C = 5/ < ADi > which is 
< C >= 2.76 ± 0.2 for the present data (galaxies with 22.25 <R< 24.0 are used 
for the shear analysis, see below). This correction does not take into account the 
presence of stars which will further dilute the signal, however, this will be minor 
for the magnitude range considered. 



4.4. Cluster Galaxies 

Including cluster galaxies in the shear measurment will reduce the measured 
value of the shear (assuming cluster galaxies are randomly aligned) by an amount 
equal to the contamination fraction. This effect is likely to be a function of 
radius since we expect the cluster galaxies to be centrally concentrated within 
the cluster. For the shear analysis in this paper we use galaxies in the range 
22.25 <R< 24.0. In order to get an estimate of the number of cluster galaxies 
in this brightness range we assume that the contamination is zero at the edge of 
our image. We then simply measure the density of galaxies in our sample as a 
function of radius and fit a straight line. This analysis reveals that at a radius of 
27.5" about 22% of the galaxies in our field sample are probably cluster galaxies 
falling to zero (by construction) at 275". We adopt this correction for the rest of 
the paper. If there are a significant number of cluster galaxies at the edge of the 
image then this correction will be too small and we will be underestimating the 
shear due to the cluster. 
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4.5. Measurements 

The value of the mean distortion, < D{f) >, for r equal to the position of 
the central dominant galaxy (CDG) in MS 1224+2007 for 874 galaxies in the 
combined R image having 22.25 < R< 24.0 in the radial range 27.5" < r < 275" 
is < D raw >= 0.034 ± 0.006; the signal-to-noise is 5.5. After correction for 
blurring by the PSF and shear polarizability (see §(4.3|) the value is < D cor >= 



0.094 ± 0.017 For comparison, the value of < D(r) > for r at the CDG location 
for the 87 PSF stars in the corrected image is 0.0009 ± 0.002, consistent with 
zero and less than 3% of < D raw >. Given that the PSF-induced distortion on 
the resolved galaxies is actually smaller than this, the residual PSF anisotropy 
has little effect on our cluster mass estimate. 

Figure ^ shows < D(f) > (corrected, for seeing, shear polarizability and 
contamination by cluster galaxies) vs. projected radius for r at the CDG position. 
The shear profile is very flat, the best fit singular pseudo-isothermal sphere (SPIS) 
shown in Figure |2] is ruled out at the 95% confidence level (x 2 = 17.07, nine 
degrees of freedom). If we restrict ourselves to nonsingular pseudo-isothermal 
spheres (NPIS) of form: 



m) ~ [l + (R/Rc)*]W (4) 

then the best fit model has a core R c = 45±i" = 135±21^ _1 k P c (x 2 = 10-43, 
nine degrees of freedom), and a central surface density k(0) = 0.30 ± 0.055. The 
95% confidence limits on r c are 60 - 400 h^ 1 kpc. 

In order to convert the density profile into physical mass units we must know 
the value of X cr jt which requires knowledge of the redshift distribution of the 
background galaxies. There is no current redshift survey complete to R = 24.0 
(actually R = 23.9 because of reddening). We can get an upper limit on E cri t by 
using the redshifts from the Canada France Redshift Survey (CFRS) ( |Lilly et aT. 



1995| ). We derive R-band magnitudes by simple averaging of their V and I pho- 
tometry which should be accurate to 0.1-0.2 mag flFukugita et al. 1995| ). Integrat- 
ing £ cr it over the redshift distribution for 22.25 < R < 24.0 yields (E cr jt) = 1.50/j 
(Qq = 1.0) and a mean redshift (z) = 0.69. The value of Yi crit will be an overesti- 
mate since the CFRS survey is incomplete for much of this magnitude range. We 
perform a similar test with the deeper but smaller Hawaii Deep Field Redshift 
Survey (interpolating R = 1/3B + 2/31) and obtain {T> crit ) = 1.49/t and a mean 
redshift (z) = 0.69, although once again this survey is not complete over the en- 
tire magnitude range. Another approach which should be somewhat less sensitive 
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Fig. 2. — A plot of the distortion, < D >, vs projected radius in radial bins containing 87 
galaxies each (outermost bin has 91) with 27.5" < r < 275" and 22.25 < R < 24.0, centered 
on the dominant cluster galaxy. The points are from the data (la error bars) and the lines 
are the expected distortion for two mass models. The dashed line is a SPIS while the solid 
line is a NPIS with a core radius of r c = 45" = 135/i -1 kpc. A recovery factor of C = 2.76 
(see text) has been applied along with a correction for contamination by cluster galaxies. 
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to incompleteness at the faint end is to simply calculate the value of E crii for an 
object at z = 0.69 (which is the mean of both the redshift surveys). This yields 
E cr j t = 1.35/i. Finally we can use theoretical models of galaxy formation and 
evolution to extrapolate the known redshift distribution to fainter magnitudes. 
Adopting the model of |Gronwall (1996)1 we find (z) = 0.73 and (S crit ) = 1.44/i 
for the relevant reddening-corrected magnitude range. Given the numbers above 
we adopt S crit = 1.44/i with an uncertainty of around 10%. 



5. Cluster Mass 

5.1. Mass reconstruction 

In the weak lensing regime, where k « 1, the formula for the surface mass 
density is ( [Kaiser fc Squires 1993| ): 



1 ^ W(Ax,Ay,s)D t {r) 
nvr ^ Axi 2 + Ay? 

where N is the number of galaxies and n is the number density of galaxies. Eqn. 
|5] assumes that the galaxies are intrinsically (in the absence of lensing) randomly 
aligned. W is a smoothing kernel which is required to prevent infinite formal 
error. In this paper we use a smoothing kernel of the form ( [Seitz fc Schncidci 



19951) 



W( Ax, Ay,s) = l-(l+ AX + 2 AZ/ j e -i^y^ , ( 6 ) 

where V is referred to as the "smoothing scale". The variance in the dimension- 
less surface mass density is: 

<- 2 > = (7) 

2-d mass maps of MS1224+2007 are shown in Figure || and are discussed further 
in Q. 

Because of the smoothing kernel, plus biases introduced by edge effects in 
the images, Eqn || is mainly useful for determining the 2-d shapes of mass dis- 
tributions. A less biased way of obtaining mass estimates as well as azimuthally 
averaged density profiles is: 



y»2 



.. _, .. .„ x- [l-^lfl-v/i-ACT 

K[r < Ti) - K{ri <r <r ) = — ' 

10 ri < r < r . o 



N t0r .t<r A(r)(Ax,2 + Ay?) 



-10- 




Fig. 3. — Mass maps derived using Eqn. |5| with smoothing scale s = 60" (top) and s = 30" 
(bottom). A total of 1201 galaxies with 22.25 < R < 24.0 are used in these reconstructions. 
The contours are spaced in la intervals. The peak of the mass distribution is consistent with 
the position of the central dominant galaxy. North is up and East is to the left. The field is 
9.6' on a side. 
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where N io is the number of galaxies between rj and r Q . This is similar to the form 
employed by [Fahlman et al. (1994)| but is valid when k is not vanishingly small. 
Since k appears on the right hand side of the equation, an iterative approach must 
be used to obtain the density profile. Radial mass profiles for MS1224+2007 are 
shown in Figure |] and are discussed further in § [5.3| . 

Galaxy distortion is insensitive to flat sheets of mass. Consequently, all mass 
measurements described in this paper are uncertain by an unknown additive 
constant. If there is a substantial flat component to the mass distribution our 
mass estimates will be lower limits. 



5.2. 2-d Mass and Galaxy Maps 

The 2-d, KS mass map for MS1224+2007 is shown in Figure |3| superposed 
on the R-band image of the field. This reconstruction uses 1201 galaxies with 
22.25 < R< 24.0. Before discussing the mass distribution we check for potential 
biases introducecd by residual PSF anisotropies by making a map with the 87 
stars from the corrected image. The highest value in this map is 8% of peak value 
in the cluster map (equal to the l.Oer noise level of the mass map) and therefore 
all but the lowest couple of contour levels should be unaffected by residual PSF 
anisotropy. For comparison a map made from the same stars on the uncorrected 
image had a peak value five times higher. 

The mass distribution is peaked very near the position of the CDG and ex- 
hibits an extension towards the south. In the higher resolution map this extension 
appears as a distinct subcluster although its proximity to the edge of the frame 
makes this conclusion uncertain. An additional subcluster appears to the west 



of the CDG. This mass maps differs somewhat from the map of [Fahlman et al. 



(1994)| . The peak of their mass distribution is centered east of the CDG and 
appears to be extended to the North. 

In Figures |^ and |^ we show maps of the galaxy number density and flux- 
weighted number density, respectively for two R-band magnitude ranges (smoothed 
to match the high- resolution mass map). The bright number density map is 
peaked slightly to the east of the CDG while the faint map seems to have a 
concentration to the southwest of the CDG and a marginal signal in the south. 
Both show evidence for significant galaxy clustering in the northern portion of the 
field. The bright flux-weighted map appears to be bimodal with two similar sized 
luminosity concentrations. One is centered very close to the CDG (not surpising 
given the brightness of the CDG) and one is to the west. The faint flux-weighted 
map looks a lot like the faint number density map with a concentration to the 
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Fig. 4. — The upper panel is the radial mass density profile (points) for MS1224+2007 from 
Eqn. |8] assuming k = (r = 275" = 790/i~ 1 kpc). It is centered on the central dominant 
galaxy. The points are the data for the cluster derived from 874 galaxies having 22.25 < 
R< 24.0. The lower panel is the radial density profile from Eqn. [S] using k derived from 
the best fit to the shear shown in Fig. |2|. The value of r Q is 275". A recovery factor of 
C = 2.76 (see text) has been applied. The solid lines are the best NPIS model fits to the 
shear measurments shown in Figure 0. 
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southwest of the CDG and a significant feature to the north and a less significant 
feature to the south. 

It is difficult to draw strong conclusions from these maps without redshift 
information. However, the bimodal central feature seen in both of the bright 
galaxy maps and in the mass map seems to indicate that there is significant 
subclustering within the cluster. Furthermore, there is strong evidence for a 
large concentration of galaxies to the north. Because this concentration is seen 
more prominently in the faint galaxies and is not seen (or marginally seen) in the 
massmap it is likely that it is at a higher redshift than z = 0.325. 



5.3. Mass and Luminosity Profiles 

In the upper panel of Figure f| we show the azimuthally averaged surface 
mass density profile centered on the CDG as derived from Equation [S] using 874 
galaxies having 22.25 <R< 24.0. We have assumed that k = on the right hand 
side of the equation. The bottom panel shows the density profile corrected using 
the best fit values of k from the fits of a NPIS to the shear (see §^). The data 
have been corrected for seeing, shear polarizability and contamination by cluster 
galaxies. 

Figure [7| shows the mass density profile overplotted with the rest frame R-band 
surface brightness profile. The surface brightness profile includes light in galaxies 
with R > 17.28 (the magnitude of CDG) but does not include possible diffuse 
light unassociated with galaxies. It is plotted as a density contrast for comparison 
with the mass density profile. Aside from the innermost regions where the light 
is completely dominated by the the CDG, it appears that mass traces light. 
We have applied a K correction to the R-band photometry assuming the cluster 
light is dominated by early type galaxies of K cor = 0.4 ( Poggianti 1997|) and an 
extinction correction of Ar = 0.11 ([Schlegel et al. 1998| ). The mass and light 
profiles are flatter than has been seen for other clusters. Of course, subclustering 
will cause the azimuthally averaged mass and light profiles to flatten out. 

Figure |8| shows the cumulative mass and luminosity profiles and Figure ^| 
shows M(r)/L/j(r) as a function of radius. Within 100 h^ 1 kpc of the CDG the 
M/L is typical of cluster galaxies (225 ± 50 (M/L R ) Q ). However, beyond the 
region where the CDG light dominates, the M/L increases rapidly and levels off. 
Ignoring the first and last points, M/L R = 640 ± 150 (M/L B ) (M RQ = 4.31) 
which is an unusually high value. Assuming the cluster is dominated by early- 
type galaxies, the evolutionary correction to z = is -0.363 mag QPoggianti 1997] ) 
which yields M/L R = 890±200 (M/Lr) q . The typical color of a nearby early-type 
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Fig. 5. — Galaxy number density for galaxies in the range 17.28 < R < 22.25 (upper panel) 
and 22.25 < R < 23.75 (lower panel). Both have been smoothed with a Gaussian having 
a scale of 30". The contours are spaced in la steps which corresponds to 1.2 and 1.3 per 
square arcmin for the upper and lower panels, respectively 



Fig. 6. — Galaxy luminosity weighted number density for galaxies in the range 17.28 < R 
< 22.25 (upper panel) and 22.25 < R < 23.75 (lower panel). Both have been smoothed with 
a Gaussian having a scale of 30". The contours are spaced in la steps which corresponds 
to 21.5 mag per square arcsecond (R-band reddening corrected) and 24.5 mag per square 
arcsecond for the upper and lower panels, respectively 
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Fig. 7. — Plot of projected cluster mass density (points) and projected rest-frame R-band 
luminosity density (solid line) in galaxies. Both are plotted as density contrasts. The dotted 
line is the best pseudo-isothermal plus core model fit to the shear measurments shown in 
Figure |2|. 
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galaxy is (B-R) = 1.57 flFukugita et al. 1995|) which gives an evolution-corrected 



M/L B = 1300 ± 300 (M/L B ) (M BQ = 5.48). 

5.4. Comparison with Previous Lensing Measurement 

In a previous study of MS1224+2007 QFahlman et al. (1994)0 a projected 



mass estimate is given for r < 2.76' of M(r < 2.76') = 3.5 x lO 14 /^ 1 M Q (no 
uncertainty is specified but it is at least 20%). This value is obtained without 
correction for the weak lensing approximation, cluster galaxies, or mass in the 
control annulus which means it is an underestimate. The mass we find with 
the same assumptions is M= 3.6 ± 0.9 x 10 14 /i _1 M Q , in excellent agreement 
with the previous measurements. There is one caveat to this comparison in that 
different values of r Q (see Equation |]) were used for the two estimates; r Q = 400" 
was used in the previous study and r Q = 275" is used in this study. If we 
correct our value to r = 400" by extrapolating the NPIS model it increases 
the mass by approximately 20%, still within la of the previous value. If we 
correct the raw value for the weak lensing approximation it reduces the mass 
by only 3%, and correcting for all the previously mentioned systematics yields 
M(r < 2.76') = 6.3 ± 1.5 x lO 14 ^" 1 M . 

We are also consistent with the value of M/L = 800 M Q /L Q (no bandpass 
specified) determined by |Fahlman et al. (1994)| even though the cluster luminos- 
ity is calculated quite differently in this paper. In this paper the cluster M/L was 
determined by measuring the mass and luminosities as density contrasts centered 
on the CDG (see § |5.3| ). This means that we include light emanating from clus- 
tered galaxies at redshifts other than that of the main cluster provided that they 
are located near the cluster center (a uniform background will be subtracted us- 
ing this technique). In [Fahlman et al. (1994)| the luminosity estimate was based 



on 30 spectroscopically confirmed cluster members from |Carlberg et al. 1994| out 



of a sample of 75 measured redshifts to r=22.0 in a 7' x 9' field centered on the 
CDG. 



5.5. Comparison With Virial Mass Estimate 

In addition to the previous lensing mass estimate, there is a virial mass esti- 
mate based on redshifts for thirty galaxies with radii 0" < 275". The velocity 
dispersion is 802 ± 90 km s _1 flCarlberg et al. 1996| ). The velocity dispersion 



implied by our mass profile is dependent on the nature of the phase space distri- 
bution function (DF) and the extent of the cluster. If we assume an isotropic DF 
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Fig. 8. — Cumulative mass (points) and rest-frame R-band luminosity (short dashes) pro- 
files for MS1224+2007 based on vrr 2 [I(< r) - L(r -> r )\ and vrr 2 [S(< r) - S(r -> r„)\, 
respectively (see Figure |7]). These are lower limits on the true cumulative profiles. The solid 
line is the best-fit NPIS model. The dotted line shows the true cumulative mass profile if the 
the NPIS model is a good description of the data. The long-dashed line is the cumulative 
luminosity profile assuming zero background. This is an overestimate since field galaxies 
must have some contribution to the R-band light. 
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Fig. 9. — Cumulative mass-to-light ratio as a function of radius. Excluding the innnermost 
and outermost points one gets M/L R = 640 ± 150 (M/Lr) . 
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and an infinite cluster then the implied dispersion is about 1380 km s _1 averaged 
over the positions of the galaxies with measured redshifts implying a mass over 
three times higher than the virial mass estimate. If we truncate the mass profile 
at 3 h^ 1 Mpc then the projected velocity dispersion falls to 1290 km s _1 , still 
much higher than the measured value. 

There is evidence from both the mass map and galaxy maps for subclustering 
within the cluster at z=0.325 and possibly at other redshifts along the line-of- 
site. The redshift study of |Carlberg et al. 1994| identified two groups of galaxies 
at z — 0.22 and z = 0.412 for which they have estimated velocity dispersions 
of 500 km s _1 and 400 km s _1 , respectively. Subclustering within the cluster 
complicates the virial analysis while clustering along the line-of-site complicates 
the lensing analysis. For example, if the main cluster actually consists of two large 
subclusters falling together in the transverse direction then the virial analysis 
will quite likely underestimate the mass. If there are groups in the foreground or 
background then the lensing mass for the z = 0.32 cluster will be overestimated, 
although the effects on the mass-to-light estimates will depend on the redshifts 
of the groups relative to the main cluster. Subclustering at the cluster redshift 
will not effect the lensing M/L estimate. 

The current number of measured redshifts is insufficient to attempt to quan- 
titatively disentangle the contributions from the various groups and clusters. A 
full discussion will have to wait until we have completed a photometric redshift 
survey of a large region centered on the cluster. 



6. Conclusions and Future Work 

In this paper we study the mass and light distributions in the z = 0.325 
cluster MS1224+2007 out to a projected radius of 800 hr 1 kpc by measuring the 
gravitationally-induced distortions of background galaxies. We detect a shear 
signal in the background galaxies in the radial range 27.5" < r < 275" significant 
at the 5.5cr level. The resultant mass map (smoothed on 60" angular scales) 
exhibits an 8a peak centered on the dominant cluster galaxy. A higher resolution 
(30") mass map reveals evidence for substructure which is even more strongly seen 
in the distribution of galaxies. The lensing and redshift data combined indicate 
that there is substructure at the cluster redshift and galaxy clustering in both 
the foreground and background of the main cluster. 

Assuming all the detected shear is due to mass at z — 0.325 we find that, 
except in the very central regions where the light from the CDG dominates, the 
azimuthally averaged mass and light profiles follow one another with a reddening 
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and k-corrected mass-to-light ratio of M/Lr = 640 ± 150 (M/Lr) q . The profiles 
are quite flat compared to other clusters, disagreeing with a singular pseudo- 
isothermal sphere at the 95% confidence level. The best fit non- singular pseudo- 
isothermal sphere has a core radius of r c = 135l45/i _1 kpc. The flat profile is 
probably a consequence of substructure. 

Our mass and M/L estimates are consistent with the previous weak lensing 
result of [Fahlman et al. (1994) . 



If we assume that the cluster is described by a an pseudo-isothermal sphere 
with a non-singular core and has infinite extent then the velocity dispersion 
implied by the lensing mass is almost 1400 km s _1 . Truncating the profile at 
smaller radius reduces this dispersion; however, unless one truncates at very 
small radius the lensing derived value remains much higher than the measured 
dispersion of 802 km s _1 ( |Carlberg et al. 1996| ). A partial explanation might be 
subclustering for which there is strong evidence. Infall of two or more subclusters 
(approximately) perpendicular to the line-of-site could result in the virial mass 
estimate substantially underestimating the cluster mass. 

This cluster remains an anomalous object. It appears to be the highest M/L 
cluster known and therefore is a candidate for a dark mass lacking baryonic 
matter. This interpretation severely tests cosmological models which are unable 
to produce such variations in baryonic fraction on cluster scales. The conclusion is 
weakened by a lack of information regarding foreground and background clusters 
which would result in an overestimate of the cluster mass and M/L. The obvious 
next step in the study of this interesting and controversial object would be to 
increase the number of galaxies with redshift measurements. This will allow us to 
look for galaxy clustering in redshift space and disentangle the contributions of 
foreground and background mass concentrations. We can then derive definitive 
mass and mass-to-light estimates for the primary cluster and whatever large mass 
concentrations exist along the line-of-sight. The most efficient way to achieve 
these goals is with photometric redshifts. 
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